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ABSTRACT 

A  general  study  is  given  of  plans  and  cylindrical  wave  propagation 
of  combined  stresses  in  an  elastic -plastic  medium.  The  coefficients  of . 
the  governing  differential  equations,  when  written  in  matrix  notation, 
are  symmetric  matrices  and  can  be  divided  into  sub -matrices  each  of  which 
has  a  special  form.  The  relations  between  the  stresses  on  both  sides  of 
an  elastic -plastic  boundary  are  derived.  Also  presented  are  the  restric¬ 
tions  on  the  speed  of  an  elastic -plastic  boundary. 
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1.  Introduction. 


The  equations  of  motion  for  a  continuum  body  are 


bv 
St"  ' 


(i,j=l,2,3) 


(1) 


where  cr_^  is  the  stress,  the  velocity  and  p  is  the  mass  density 
of  the  body,  and  summation  is  implied  by  repeated  indices.  The  relation 
between  the  strain  and  the  velocity  v^  is 


(2) 


while  the  stress-strain  relation  for  an  elastic,  isotropic  work-hardening 
material  is  (see  [1]) 


E  is  Young's  modulus,  v  is  Poisson's  ratio,  k  is  the  yield  stress 
and  the  yield  condition  can  be  written  as 

f^j)  =  k2  (4) 

G(k)  in  Eq.  (5)  is  a  given  function  of  k  which  characterizes  the 
work-hardening  property.  Equations  (l)  ~  (4)  give  a  complete  descrip¬ 
tion  of  wave  propagation  in  three-dimensional  elastic -plastic  media. 

In  this  paper,  we  will  restrict  our  attention  to  special  cases 
in  which  the  governing  equations  depend  on  only  one  space  variable. 
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Plane  wave  propagation  and  cylindrical  wave  propagation  are  such  special 
cases.  The  most  general  plane  wave  propagation  is  the  one  in  which  v^, 
Vg,  and  v^  are  functions  of  x^  and  t  only.  Using  x  instead  of  x^ 
for  simplicity,  c^,  xQ)  x ^  for  or  respectively,  Eq.  (l) 

gives 


T3,x  =  p  v3,t 


where  the  subscripts  x  and  t  denote  partial  differentiation  with 
respect  to  these  variables.  If  we  use  von  Mises  yield  condition: 


(9) 


v..  =  -  t,  .  +  2t,GQ. 

5*x  U  3,t  3 


where  s^,  s2  and  stand  for  s.^,  s 2g,  and  respectively.  M 
is  the  shear  modulus  and 


8  '  +  2E2°2,t  +  2V2,t  +  2Vj,t 


Equations  (5  ~  ll)  give  a  complete  description  of  a  general 
plane  wave  propagation  in  an  elastic -plastic  medium. 

When  v2  s  o,  and  hence  Tg  s  0,  Eqs.  (5)*  (7)/  (8),  (9),  and  (ll) 
reduce  to  the  case  of  pressure-shear  wave  propagation  considered  in 
[2,3,4] .  If  v5,  v^  are  the  only  non-zero  velocity  components,  then 
and  are  the  only  non-zero  stresses  and  Eqs.  (6),  (7)>  (10) 
and  (ll)  reduce  to  the  case  of  two-shear  waves  studied  in  [3,5*6].  The 
equations  derived  in  [7]  for  combined  longitudinal  and  torsinal  waves 
in  a  thin-walled  tube  can  also  be  reduced  from  Eqs.  (5)*  (6),  (8)  and 
(10)  by  letting  0. 

If  v  ,  v„,  v  are  the  velocity  components  of  a  particle  in 
r  9  z 

cylindrical  coordinates  (r,0,z),  the  most  general  cylindrical  wave 

propagation  is  the  one  in  which  v  ,  v^  and  v^  are  functions  of  r 

and  t  only.  Then,  the  only  non-zero  strains  are  er,  e0,  70*  7Z 

where  7  stand  for  e  e  respectively.  Consequently  cr , 

'  v  *  z  tv  rz  * 

ct  ,  o  ,  t  ,  t  are  the  non-zero  stresses  where  t  --  cr  ,  t  =  cr. 

0  z  6  z  v  tv  z  tz 

Now,  instead  of  writing  the  governing  equations  for  cylindrical  wave 
propagation  in  the  form  shown  in  Eqs.  (5  ~  ll)*  we  will  use  matrix 
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notation  and  write  the  governing  equations,  Eqs.  (l~3)  in  a  matrix 
differential  equation  which  can  he  applied  to  both  plane  and  cylin¬ 
drical  wave  propagation. 

2.  The  Matrix  Differential  Equation. 

For  wave  propagation  which  involves  only  one  space  variable  x 
or  r,  the  equations  of  motion  (l)  can  be  written  in  matrix  notation 
as 


M  cr  +  b,  *  p  v.  (15) 

- x  ~1  ~t 

where  cr  and  v  are  column  vectors  whose  elements  are  stress  and 

♦>/  /v 

velocity  components  respectively,  b^  is  also  a  column  vector  whose 
elements  are  functions  of  stress  or  and  space  variable  x  only,  (in 
cylindrical  waves,  x  becomes  r).  If  cr  has  m  elements  and  v 
has  n  elements,  then  M  is  an  n  X  m  matrix  whose  elements  are 
constants.  The  continuity  condition  Eq.  (2)  can  be  written  as 


6 ,  =  N  v 
~t  ~  ~x 


where  e  is  a  column  vector  with  strains  as  its  elements  while  b0 
is  a  column  vector  whose  elements  are  functions  of  velocity  v  and 
space  variable  x  only.  N  is  an  m  X  n  matrix  whose  elements  are 
constants.  Finally,  the  stress-strain  relation,  Eq.  (2)  can  be 
written  as 


e.  =  S  cr 


(15) 


where  S  is  an  m  X  m  square  matrix .  In  the  plastic  region,  S 
can  be  written  more  precisely  as 
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SP=  Se  +  G(k)(yf)(vf)T  (l6) 

0 

S  is  also  an  m  X  m  square  matrix  whose  elements  are  functions  of 
elastic  constants  only.  Thus 


(17) 


gives  the  elastic  stress-strain  relation.  Vf  is  the  gradient  of 
f(g)  with  respect  to  the  components  of  a.  Hence,  by  (4) 


(Vf)T  S,.  =  2kkt 


(18) 


Now,  by  eliminating  between  Eqs.  ( 14)  and  (15),  we  can  write 
Eqs.  (13)~(15)  in  one  matrix  equation: 


where 


(19) 


(20) 


governing  equations  for  general  plane  wave  propagation  derived  in 
Eqs.  (5~ll)  as  well  as  the  equations  for  other  plane  wave  propagation 
reduced  from  Eqs.  (5~ll)  can  be  written  in  the  form  of  Eq.  (19) • 


Moreover,  A  and  B  are  both  symmetric  in  all  cases.  In  particular, 
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T 

N  =  M  and  the  elements  of  N  and  M  are  either  one  or  zero. 

The  same  is  true  for  all  cylindrical  waves.  For  the  most  general 
cylindrical  wave  propagation,  it  can  be  shown  that,  using  von  Mises 
yield  condition: 


V 

r 

(^-OoVr 

{ 

'1  0  0  0  0 

11 

i>l 

ye 

V 

Z 

'  *1- 

2Vr" 

Vr 

,  M  = 

0  0  0  1  0 

.0  0  0  0  1 

N  = 


(21) 


l/E 

-v/E 

-v/E 

0 

0 

-v/E 

l/E 

-v/E 

0 

0 

-v/E 

-v/E 

Vs 

0 

0 

II 

0 

0 

0 

Vn 

0 

0 

0 

0 

0 

V(i 

Cristescu  [8]  has  studied  an  axilly  symmetric  wave  propagation 

in  which  vr  and  v0  are  the  only  non-vanishing  velocities .  The 

governing  equations  can  be  obtained  from  Eq.  (21)  by  letting 

v  =  t  =  0.  The  particular  cases  of  two-shear  cylindrical  waves 
z  z 

and  pressure-shear  cylindrical  waves  have  also  been  derived  earlier 
in  [5]  but  the  resulting  coefficient  matrices  were  not  symmetric. 

The  formulation  presented  here  yields  symmetric  coefficient  matrices 
for  all  cylindrical  waves. 
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The  analyses  in  the  rest  of  this  paper  will  be  based  on  the 
matrix  differential  equation  ( 39) •  Although  for  all  plane  waves 
and  cylindrical  waves  the  matrices  A  and  B  are  symmetric,  the 
analyses  presented  in  the  following  do  not  require  the  symmetry 
property  of  A  and  B.  The  yield  condition  f  of  Eq.  (4)  is  not 
restricted  to  von  Mises  yield  condition. 
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Now,  by  the  theory  of  determinants,  it  is  known  that 


llp-i +  a  h  ,  p  ,  p,,-..>pjl=lip||  +  a||h,  d^,  p,,...,pj|  ,  (26) 

ao  o'*  ^  o"-  aoJ* 

||h,  ah,  p,,  pk,  •  *  - , P J|  =  0  (27) 

By  repeatedly  applying  Eqs.  (26)  and  (27)  to  the  right  hand  side  of 
Eq.  (25),  we  obtain: 

112  +  a£  fi1®  =  ®Ei'  Ey'Erll  +  a  8i*£-  Ea’  Ey-'Er11 

(28) 

+  a  gjp,,  h,  p,^  •  •  •>  pjl  +  +  a  gjp,,  . . .  ,fa|| 

With  this,  it  is  not  difficult  to  see  that  Eq.  (28)  can  be  written  in 
the  form  of  Eq.  (22).  This  completes  the  proof. 

For  the  particular  case  in  which  P  is  a  unit  matrix,  Eq.  (22) 
reduces  to 


I  +  a  h  gT||  =  1  +  a  hTg 


5.  The  Characteristic  Equation. 

The  characteristics  c  of  Eq.  (19)  are  ihe  roots  >f  the  equation 
(see  [9]): 

|| cA  -  B||  *  0  .  (50) 


Since 


cA-B  = 


pci 

ao 

M 

N 

pci 

1  AO 

o' 

AO 

*  I 

-i-  M 
pc  ~ 

N 

o' 

cS 

AO 

N 

AO 

■ 

D  ’ 

AO 

> 

? 

:  0 

1  ao 

\ 

V. 

—  I 
pc 
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where 


D  =  pc2s-Nfc$,  (31) 

we  have 

l|cA-B||  =  (Pc)n'my!  •  (32) 

Thus,  instead  of  expanding  the  determinant  ||cA-b||  which  is  of  order 
m+n,  we  can  expand  the  determinant  ||g!l  which  is  of  order  m. 

If  we  define 

De  =  pc2ge‘NM,  (33) 

then  by  Eos >  (l6)  and  (31),  we  have 

DP  =  De  +  pc2G(k)(Vf)(yf)T  •  (34) 

Using  the  identity  derived  in  the  previous  section,  we  obtain 

lljfll  =  llse||  +  Pc2G(k)fef)T(~e)*(Vf)  .  (35) 

Equation  (35)  can  be  used  to  study  the  relative  positions  of  the  roots 

of  IIbPII  =  0  and  ||jDe(|  =  0.  (see  [4]). 


4.  The  Elastic -Plastic  boundary. 


If  c  is  the  speed  of  an  elastic -plastic  boundary,  then 


dw 

—  =  w  c  +  w , 
dt  ~x  ~t 


(56) 


is  the  total  derivative  of  w  along  the  boundary.  Elimination  of  w 
between  Eqs.  (19)  and  (36)  yields 
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(57) 


(°A-B)wt  =  cb-g^  . 


Since  w  is  continuous  across  an  elastic -plastic  boundary,  Eq.  (57) 
gi\  es 

(cAP-  S)w^  =  (eAe  -B)w^  (5 

where  the  superscripts  e  and  p  denote  the  values  in  elastic  and 
plastic  regions  respectively.  With  Eq.  (20),  Eq.  (58)  is  equivalent 


p  c  vp  +  M>r?  =  pcv?  +  M  of 
~~t  *t  ~~t 


N  y?  +  c  SP£?  =  N  vf  +  c  Se£^  • 


It  is  more  convenient  to  write  Eqs.  (59)  and.  (4o)  in  the  following  forms, 
using  Eqs-  (l6)  and  (l8): 


pc(v^  -vjl)  +  M(o^  -o£)  =  0 


N(^-yP)  +  =  2cG(k)k  kP(Vf) 


v^t  H 


(42a) 


N(v^-yP)  +  cgP(^-o£)  =  2cG(k)k  k*(Vf) 


'-H  ~t 


(42b) 


Equations  (4l)  and  (42a)  yield 


Se(S^-4)  =  2pc2G(k)k  kP(Vf) 


(45a) 


while  Eqs.  (4l)  and  (42b)  yield 


BP(aJ"S^)  =  2pc2G(k)k  k®(Vf) 


(45b) 
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Equation  (43a)  or  (43b)  gives  a  relation  between  cr^  and  on 
both  sides  of  an  elastic-plastic  boundary-  For  the  case  of  combined 
longitudinal  and  torsional  stresses  in  a  tube., 


D6  = 


r  2 
pc 

E 


-  1 


pc 


-  1’ 


and  Eq-  (43a)  reduces  to  the  result  obtained  in  [10]. 

Equation  (43a)  can  be  solved  for  (o^-£^)  by  pre -multiplying  both 


sides  of  equation  by  (De)  ^  and  _aking  use  of  Eq.  (23).  Hence 


II  On  T.'^  '  ^ 


M 


(44) 


m 

Pre -multi plying  once  more  both  sides  of  Eq.  (44)  by  (£f )  and  making 
use  of  Eq.  (l8),  we  obtain 


**  (Vf)T(Der(W) 


l|De|| 


(45) 


Since  the  right  hand  side  of  the  equation  is  a  scalar, 

(Vf)T(De)*T(Vf)  =  (Vf)T(De)#(Vf), 
and  by  virtue  of  Eq.  (35)  we  have  finally 


fc*  ll2PH 


(46) 
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A  simplest  particular  case  of  Eq.  (46)  is  the  case  of  longitudinal  wave 
propagation  in  a  thin  rod  studied  in  [ 13 3 •  In  the  case  of  combined 
longitudinal  and  torsional  stresses  in  a  thin-walled  tube, 


||De||  =  (%  '  1)(4  " 

c  c„ 


iidpii  =  A  - 1 )A  - 1) 


where  c2  =  E/p  ,  c2  =  p/p  ,  and  c  ,  c  are  the  roots  of  ||dp||  =  0, 

O  u  IS 

Eq.  (46)  reduces  to  the  result  obtained  in  [ 10] - 

If  kP  =  0,  then  by  Eq.  (43a)  either  or  ||De||  =  0. 

Similarly,  if  k®  =  0,  then  by  Eq.  (43b)  either  a£  =  or  ||DP||  =  0. 

Consequently,  if  k®  =  kp  =  0,  Moreover,  v^  =  yp  by  Eq. 

(4l).  Thus  when  k®  =  kP  =  0,  *  wj.  To  obtain  relations  between 

ae  and  qF.  on  both  sides  of  an  elastic-plastic  boundary  when 
~tt  ~tt 

kf  =  kp  =  0,  we  differentiate  both  sides  of  Eqs-  (43)  and  noticing 

"fc  u 

that  £  =  4  ,  kt  =  kt  =  0, 


De(4-aPt)  =  2P02G(k)kkPt(Vf) 


(47a) 


dp(4  -<&.)  =  2Pc2G(k)k  kjt(52f)  • 


(47b) 


Similarly,  if  we  differentiate  both  sides  of  Eq.  (45),  we  obtain, 
after  using  Eq.  (35), 

4  hfl 
4  ‘  tel ' 
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Equations  (47)  and  (48)  are  identical  to  Eqs.  (43)  and  (46)  respectively 
with  the  exception  that  the  order  of  derivatives  is  changed. 

From  the  above  derivations  we  can  generalize  the  restilts  and  state 
as  in  the  following: 

If  Snke/Stn  =  Snk^/ Stn  =  0  for  n=l,2, . .  .,<2-1  but  cPke/ Sta 
and  are  not  both  zero,  we  have 


or 


\n  e  \n  p 

o  w  o  wr 


De( 


St” 

b%e 

Sta 

St0' 

■>a  e 
o  a 

1?' 

Sta 

,  n— 1, 2,  •  * ♦  ,Q:-1 


)  =  2pc<1G(k)k  (£f) 

St 


32G(k)k  ^ 

Sta 


(49) 


(50a) 


(50b) 


and 


aV/ata  W 
Sakp/3t“  '  l|DeIJ 


(51) 


For  a  loading  wave, 


oake/Sta  J  >  0 

sV/sta  \  <  0 


if  d  is  an  odd  integer 
if  a  is  an  even  integer 


(52a) 


while  for  an  unloading  wave, 

oake/Sta  f  <  0 

sV/sta  \  >  0 


if  a  is  an  odd  integer 
if  a  is  an  even  integer 


(52b) 
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With  this,  Lq.  (51)  furnishes  a  restriction  on  the  speed  of  an  elastic  - 
plastic  boundary  as  in  the  particular  cases  considered  in  [10,11,13]. 

It  whould  be  noticed  that  Eq.  (5l)  applies  if  dek/dtn  (n=l,2, . . .,a-l) 
are  zero  on  both  sides  of  an  elastic -plastic  boundary  as  well  as  along 
the  boundary.  If  dnk/dtn  (n=l,2, . . .,a-l)  are  zero  on  both  sides  of 
the  boundary  only  at  the  point  concerned  but  not  at  other  points  on  the 
boundary,  the  right  hand  side  of  Eq.  (51)  should  be  modified.  An  example 
of  this  modification  for  the  particular  problem  of  wave  propagation  in 
a  rod  was  given  in  [12]. 
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